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The hyperfine induced 2s 2p ^Pq ^ 2s^ ^5*0 transition rate in Be-like sulfur was measured by 
monitoring the decay of isotopically pure beams of 32gi2+ 33gi2+ ^^^^ ,^ heavy-ion storage 
ring. Within the 4% experimental uncertainty the experimental value of 0.096(4) s~^ agrees with 
the most recent theoretical results of Cheng et al. [Phys. Rev. A 77, 052504 (2008)] and Andersson 
et al. [Phys. Rev. A 79, 032501 (2009)]. Repeated experiments with different magnetic fields in the 
storage-ring bending magnets demonstrate that artificial quenching of the 2s 2p ^Pq state by these 
magnetic fields is negligible. 

PACS numbers: 32.70.Cs, 31.30.Gs, 34.80.Lx 



I. INTRODUCTION 

Hyperfine quenching in atoms and ions [1] is the short- 
ening of excited-state lifetimes by the interaction of the 
electron shell with the magnetic moment of the atomic 
nucleus. A particularly drastic hyperfine quenching ef- 
fect is observed in alkaline-earth-like and, in general, di- 
valent atoms and ions (with valence shell n) where the 
first excited level above the ground state is the ns np ^Pq 
state. (Fig. 1). A total electronic angular momentum 
of J = for this level makes a single-photon decay to 
the {nsY ^So ground state impossible. The nsnp ^Pq 
states are thus extremely long-lived considering the fact 
that the lifetimes associated with the two-photon ElMl 
transition can be up to several months [4, 5] . 

However, a nucleus with nonzero spin induces a mix- 
ing of the ns np ^Pq state with its neighboring ns np "^Pi 
state via the hyperfine interaction. This drastically in- 
creases the nsnp '^Pq (ns)^ ^Sq transition rate. Such 
hyperfine-induced (HFI) decay rates have been calculated 
theoretically for beryllium-like (n=2) [6-10], magnesium- 
like (n=3) [7, 11-13], and zinc-like (n=4) ions [14-16]. 
Calculations have been carried out also for divalent heav- 
ier atoms [17, 18] and singly charged ions [19], where 
the long isotope-dependent lifetimes are attractive sys- 
tems for obtaining ultraprecise optical frequency stan- 
dards [20, 21]. A comprehensive review of the field, which 
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FIG. 1. Simplified level diagram for beryllium-like S . The 
level energies [2] are given in eV on top of the horizontal lines 
which are not drawn to scale. The level lifetimes (in s) la- 
beling the excited levels on the right side of the horizontal 
lines were calculated from theoretical one-photon transition 
rates [3] that do not account for hyperfine effects. The one- 
photon transitions are labeled El (electric dipole), M2 (mag- 
netic quadrupole), and HFI (hyperfine induced). Numbers 
in square brackets denote powers of 10. In case of nonzero 
nuclear spin the hyperfine induced ''Po ^ So transition rate 
AuFi acquires a finite value. 



also covers other classes of atomic systems such as He- like 
ions, has been published recently [1]. 

So far, very few experimental measurements of 
hyperfine-induced (HFI) decay rates for ns np '^Pq states 
in divalent ions have been performed. Experiments with 
singly charged In"*" and Al"*" ions [22, 23] have used radio- 
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TABLE I. Comparison of theoretical values for the hyperfine 
induced 2s 2p ^Pq 2s^ ^Sq transition rate in Be-like Ti^®+ 
with the only published experimental result. 



Year 


^HFI (S 

Experiment 


Theory 


Deviation from 

experiment 


Ref. 


1993 




0.3556 


-37% 


[6] 


2007 


0.56(3) 






[27] 


2008 




0.6727 


20% 


[8] 


2009 




0.6774 


21% 


[9] 


2010 




0.661 


18% 


[10] 



frequency traps. The rates for the atomic-clock transi- 
tions 5s5p ^Pq -> 5s^ ^5*0 in ^^^In+ and 3.s3p ^Po 
3s^ in ^''Al+ have been measured using optical meth- 
ods with uncertainties [24] of ~ 4% and ~ 7%, respec- 
tively. For Be-like N'"'"'" the HFI rate has been extracted 
from observations of a planetary nebula [25] with an un- 
certainty of 33%. In astrophysics HFI transitions can be 
used to infer isotopic abundance ratios which, in turn, 
provide insight into stellar nucleosynthesis [26]. 

A laboratory measurement [27] of a HFI transition rate 
in a multiply charged beryllium-like ion was carried out 
with Ti-'^®+ employing the same storage-ring technique 
as in the present work. The experimental value for the 
HFI 2s 2p transition rate is 0.56(3) s'^ 

The only theoretical value [6], that was available at the 
time of the experiment, was significantly smaller (Tab. 
I). This first laboratory measurement of a HFI induced 
transition rate in a Be-like ion has triggered new theo- 
retical calculations not only for Be-like ions [8-10] but 
also for Mg-like ions [12, 13] and Zn-like ions [15]. Ta- 
ble I compares the experimental HFI 2s 2p ^Pq 2s^ ^S'o 
transition rate for Be-like 47'pji8-i- -^w^]^ presently avail- 
able theoretical values. The new theoretical results agree 
slightly better with the experimental finding than does 
the older theory. In contrast to the early value, which is 
smaller than the experimental rate by 37%, the new ones 
are larger by up to 21%. 

The main technical difference between the old and the 
new calculations is the treatment of the atomic structure. 
Be-like ions have a rather compact electron shell where 
correlation effects are particularly strong. It turns out 
that the calculated values for the 2s 2p 2s^ transition 
rates depend sensitively on the theoretical description of 
the atomic states, i.e., on the size of the configuration 
space considered and on the inclusion of relativistic and 
QED effects [9]. HFI transition rates are thus excellent 
probes of atomic correlations. 

At present, the origin for the remaining significant 
deviation by 18% from the experimental result is un- 
clear. It has been suspected that the magnetic fields in 
the storage-ring bending-magncts lead to an additional 
quenching of the 2s 2p ^Pq state. However, a recent 
theoretical investigation [28] concluded that this effect is 
negligible under the conditions of the Ti""^*"*" experiment. 

Here we present storage ring measurements of the HFI 



2s2p^Po — ^ 2s^ ^5*0 transition rate in Be-like ^^S^^+ us- 
ing the same experimental technique as before [27]. The 
aim of the present investigation is to enlarge the exper- 
imental data base and to possibly shed some light onto 
the origin of the discrepancies discussed above. In par- 
ticular, measurements have been carried out at different 
magnetic field strengths in order to experimentally inves- 
tigate the magnetic quenching effect. 

II. EXPERIMENT 

The experiment was performed using the accelerator 
and storage-ring facilities of the Max-Planck-Institut fiir 

Kernphysik (MPIK) in Heidelberg, Germany. Negatively 
charged ions from a sputter ionization source were in- 
jected into the MPIK tandem accelerator. In order to 
produce ^^S beams, isotopically enriched FeS samples 
with a relative '^^S abundance of 99.8% were used in the 
ion source. The FeS samples for the production of '^^S 
beams contained sulfur with a natural isotope distribu- 
tion (Tab. II). 

Multiply charged sulhir ions were produced by electron 
stripping in a gas stripper located on the high-voltage 
terminal of the MPIK tandem accelerator and, subse- 
quently at a higher ion energy, in a 5 /ig cm~^ carbon 
foil located behind the exit of the tandem accelerator. 
Because of the statistical nature of the electron stripping 
process, ions were produced in a range of charge states. 
The stripping also leads to internal excitation of the re- 
sulting ions. Most of the excited states decay rapidly 
to the ground state and to long-lived excited states such 
as the S^^+(2s2p ^P) states (Fig. 1) which are of par- 
ticular interest in the present work. The centroid of the 
charge state distribution depends on the ion energy which 
was chosen to be about 56 MeV at the carbon foil. At 
this energy, charge state 12+ was produced with nearly 
maximum efficiency. This charge state was selected for 
further beam transport by passing the beam through 
a dipole magnet which dispersed the various ion-beam 
components according to their mass-to-charge ratios. 

Isotopically pure pulses of S^^"*" ions were then in- 
jected into the heavy-ion storage-ring TSR. The dipole 
bending magnets of the storage ring were adjusted such 



TABLE II. Some relevant properties of 32312+ ^^^^ 33gi2+ 
ions. Listed are the mass number A, mass m [29, in atomic 
mass units], nuclear spin I, nuclear magnetic moment fj, [30, in 
units of the nuclear magneton /iat], natural abundance / [31], 
and the lifetime teimi [4, HFI quenching not considered], 
due to the 2s 2p ^Po -> 2s^ ElMl two-photon decay. The 
uncertainties of m, ji, and / are smaller than the respective 
numerical precision of the tabulated values. 
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FIG. 2. Sketch of the heavy-ion storage-ring TSR. The 8 
bending dipole magnets with 1.15 m bending radius together 
cover 13% of the closed orbit of the ions which has a circum- 
ference of 55.4 m. Both the electron target (Target) and the 
electron cooler (Cooler) were used for the present electron-ion 
recombination measurement. The locations of the recombina- 
tion detectors are indicated by the thick curved arrows in the 
top left and bottom left corners of the ring which represent 
the paths of recombined ions. The labels PGl - PG4 denote 
the approximate locations of the pressure gauges that were 
monitored by the data acquisition system. 



that the S^^+ ions were circulating on a closed orbit. 
Stored S^^+ ion currents were up to 35 /lA. The TSR 
storage ring [32] is equipped with two electron-beam ar- 
rangements, which are referred to as the "Cooler" and 
the "Target" (Fig. 2) and which can be used for electron 
cooling and for electron-ion collision studies. For the 
present measurements various modes of operation were 
employed where either the Target or the Cooler was used 
for electron cooling [33] and also as an electron target for 
electron-ion collision experiments. Details will be given 
below. 

The electron-ion collisions that occur in the merged- 
beams regions can lead to formation of S^^"*" ions. These 
recombination products are less strongly deflected in the 
TSR dipole magnets than the more highly charged pri- 
mary ions. This enables the recombined S^^^ ions to be 
easily counted with single-particle detectors which are 
positioned such that they intercept the path of the S^^+ 
ions (Fig. 2). A channeltron-based detector [34] was used 
behind the Cooler and a scintillation detector [35] was 
employed behind the Target. Because the ions moved 
with high velocities of 5-10% of the speed of light the 
detection efficiency was practically 100% for both detec- 
tors. 
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FIG. 3. Measured DR spectra of 32312+ (f^Q symbols, top 
curve, offset by -f lO^'"* cm"^ s~^) and ■'■^312+ ^Qpgjj symbols, 
bottom curve). For these measurements the Target was used 
as electron target and the Cooler was used to continuously 
cool the ion beam during recombination data taking. The 
vertical arrow marks the position of a resonance which is asso- 
ciated with primary ions in the 2s 2p ^Po state (see discussion 
in the text). 



A. Dielectronic recombination spectra 

After injection of ion pulses into the TSR and sub- 
sequent electron-cooling of the stored ion beam, the re- 
combined S^^"*" ion signal was recorded as a function of 
electron-ion collision energy for both 32gi2-i- ^^^^ 33gi2-i- 

ions using well established procedures (e.g., [35-37], and 
references therein). Further details and comprehensive 
results of the S^^"*" recombination measurements will be 
published elsewhere, since only a very narrow range of 
electron-ion collision energies is of interest for the present 
study. 

Figure 3 shows recombination spectra of 32gi2+ 
33gi2+ jQj-^g ^j^g 0.28-0.49 eV energy range. In both 
spectra two strong DR resonances appear at 0.315 and 
0.460 eV. These belong to the (2s 2p ^P)9l manifold 
of resonances associated with 2s 2p excitation of 
S"'^^"'"(2s^) ground-state ions. The much weaker resonance 
at 0.4 eV is only visible in the 32gi2-i- spectrum. It is as- 
sociated with the excitation of ions which are initially in 
the practically infinitely long-lived 2s 2p ^Pq metastable 
excited state. The resonance appears rather weak since 
only a relatively small fraction of Be-like ions were in the 
2s 2p ^Pq state. For isoelectronic Ti^^+ [38] and Fe'^^+ 
[39] ions stored in TSR, metastable fractions of 5% and 
7%, respectively, were inferred. A similar value can be as- 
sumed here. The lifetimes of the 2s 2p ^Pi and 2s 2p ^P2 
states (Fig. 1) are orders of magnitude shorter than the 
1.8 s long time interval that was used for electron cool- 
ing after each injection of new ions into the storage ring 
and before data taking was started. Therefore, it can 
be safely assumed that the fractions of primary ions in 
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the 2s 2p "^Pi and 2s 2p states were negligibly small 
during all measurements. 

For the ^^S isotope, the ^Pq resonance at 0.4 eV does 
not contribute to the DR spectrum because the 2s 2p ^Pq 
state is quenched by hyperfine interaction. There is suffi- 
cient time for its decay since the ions are stored for much 
longer than the longest predicted S^^+(2s2p ^Pq) HFI 
lifetime of about 28 s [6, 8-10]. 



B. Resonance-decay curves 

In order to obtain an experimental value for the HFI 
lifetime of the 2s 2p ^Pq state, the decay of the 0.4 eV 
resonance was monitored as a function of storage time. 
To this end only one ion pulse was injected into the 
TSR storage ring. After a short initial cooling period 
the electron-ion collision energy in the Cooler was set to 
0.4 eV and the recombination signal from the ^Pq DR 
resonance was recorded for up to ~220 s. Then, the re- 
maining ions were kicked out of the TSR and a new ion 
pulse was injected. This sequence was repeated for sev- 
eral hours in order to reduce statistical uncertainties to 
a level as low as achievable within the limited amount of 
available beam time. 

Because of a small vacuum leak in the recombina- 
tion detector at the Target, this detector could not be 
used without compromising the residual gas pressure in 
the TSR and, consequently, considerably shortening the 
beam lifetime. Thercifore, for all resonance-decay curve 
measurements the Cooler was used as an electron target. 

As mentioned above, a change of the isotope required 
the change of the sample in the ion source. Moreover, 
the accelerator settings had to be slightly retuned. The 
entire procedure for changing the isotope took 2-4 hours. 
Therefore, isotope changes were performed only every 
other day. 

Different data sets (labeled A, B, and C) were obtained 
under quite different experimental conditions. Data set A 
was measured in October 2010. Data sets B and C were 



TABLE III. Some of the experimental parameters that were 

used during the decay-curve measurements: Atomic mass 
number (A), magnetic field in the TSR-dipolcs (iJdip), ion 
velocity (TJion), relativistic 7-factor of the ions, and electron 
density (rie) in the Cooler at 0.4 eV electron- ion collision en- 
ergy. 
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FIG. 4. Large panels: Measured resonance-decay curves 
(symbols) and fits of Eq. 12 to the data (full curves). Small 
panels: Residuals from the fits. 



collected in April 2011. Some experimental parameters 
are given in Tab. Ill and further details are disciissed 
below. All measured beam-decay curves are displayed in 
Fig. 4. 



1. General properties 

In the decay curves in Fig. 4, a few general properties 
can be observed: 

(a) Substantial differences are seen between the three 
pairs of curves for the data sets A-C, respectively. This 

arises because the electron-induced count rates are super- 
imposed on backgrounds of varying sizes, which originate 
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from from charge changing coUisions of the ions with the 
residual gas in the TSR. In sets A and B, the relative 
contribution of the background to the measured signal is 
much larger than in C. Therefore, the A and B curves are 
dominated by the collision- induced ion beam loss. The 
electron-induced signal here represents a small contribu- 
tion only, though it still can be analyzed as discussed be- 
low. In contrast, the electron-induced signal dominates 
for case C. 

(b) In Fig. 4i of data set C, the electron-induced de- 
struction of primarily the mctastable S^^~^{2s2p ^Pq) 
ions leads to a rather fast beam loss during the first ~30 s. 
Later, only the collision-induced signal from the ground- 
state beam, with longer lifetime, is left. 

(c) In Fig. 4j of data set C, an even faster initial de- 
cay is clearly seen. This is caused by the additional, 
HFI radiative decay of the metastable state, leading to a 
corresponding sharp decrease of the electron-induced sig- 
nal, before the collision-induced background takes over. 
Hence, in the comparison of Figs. 4i and 4j the HFI 
quenching is very clearly visible. 

(d) In the curves for data sets A and B, this behavior 
occurs much less prominently because of the large back- 
ground. However, small fast decay components from the 
HFI decay can be seen also in Figs. 4b and 4d, comparing 
to the corresponding set for the even-A isotope. 

Using the model given in Sec. HI A, time constants 
for the various decay modes can be extracted from all 
three data sets. First, though, we separately specify the 
measuring conditions of the sets. 



2. Data set A 



Data set A (Figs. 4a and b) was taken with the ions 
stored at their injection energy of about 56 MeV. More 
precise values for the ion velocity are given in Tab. III. 
The Cooler was used as a recombination target with the 
electron-ion collision energy set to 0.4 eV. The Target was 
used for continuous electron cooling of the ion beam. 

This mode of operation provided stable reproducible 

conditions only for up to 24 hours. After that the photo- 
cathode [40] of the Target routinely had to be replaced by 
a freshly prepared one. It turned out that the emission 
characteristics of the various photocathodes used were 
significantly different, causing different Target electron 
beam profiles, and the residual gas pressures in the Tar- 
get region of the TSR changed after a new cathode was 
installed at the Target. Since the altered conditions influ- 
enced the ion-beam lifetime considerably, measurements 
before and after a cathode change could not be combined 
into one data set. Nevertheless, we were able to measure 
one set of decay curves under stable conditions (Fig. 5) 
with one photocathode for both isotopes. 
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FIG. 5. Vacuum pressures at different locations (labels PGl-4 
in Fig. 2) in the TSR storage ring recorded during the life- 
time measurements which contributed to data set A with a 
single photocathode. Time zero is October 1st, 2010 at 16:14. 
During the gap between 12 and 16 h the isotope was changed 
from to ^^S. The pressures remained nearly constant over 
the entire measurement time of data set A. It should be noted 
that the absolute pressure values bear rather large (unknown) 
uncertainties. 



3. Data set B 

In order to avoid complications due to Target cathode 
changes, the Target was switched off during the measure- 
ments of data sets B and C. The Cooler was used for both 
electron cooling and as recombination target. The cooler 
is equipped with a thermionic cathode which can eas- 
ily be operated with stable performance for several days. 
The electron energy was switched alternatingly between 
cooling energy, i.e., zero electron-ion collision energy, and 
an electron-ion collision energy of 0.4 eV. Dwell times of 
20 ms were used after each energy jump and data were 
taken only during the last 10 ms of each dwell time in- 
terval. The data that are shown in Figs. 4e and 4f are 
the data that were accumulated while the electron-ion 
collision energy had been set to 0.4 eV. 

The measurements of data set B were carried out with 
practically the same ion velocities as were used for data 
set A. However, the slope of the data set B decay curves 
(Fig. 4e and 4f) is steeper than the slope of the data set 
A curves (Fig. 4a and 4b). This change in beam lifetime 
by about 50% is largely due to differences between the 
TSR vacuum conditions of October 2010 and those of 
April 2011. 



4- Data set G 

For data set C (Figs. 4i and 4j) a higher ion energy of 

about 225 MeV was used. The higher energy was real- 
ized by operating TSR in synchrotron mode after each 



6 



injection of a new ion pulse. In order to keep the ions 
circulating during acceleration the magnetic fields were 
ramped synchronously to higher values. No more than 
50% of the ions were lost during this phase which lasted 6 
s. Ion currents after ramping amounted to up to 17 /xA. 
The magnetic field in the TSR dipole magnets was in- 
creased roughly by a factor of 2 (Tab. III). The same 
sequence of jumps in electron-ion collision energy was 
applied as for data set B. Because of the higher ion ve- 
locities this had to be realized at correspondingly higher 
laboratory electron energies. 



III. DATA ANALYSIS AND RESULTS 



A. Time evolution of state populations 



Under the assumption that the target electron density 
and the partial pressures of the residual-gas components 

do not change during the measurement, the time evo- 
lution of the populations Nm and Ng of the metastable 
2s 2p ^Pq state and of the 2s^ ^Sq ground state, respec- 
tively, can be obtained from solving the following set of 
rate equations [41]: 



where the rates Ax refer to various atomic processes that 
populate or depopulate states m and g. The rate A^. 
is the radiative HFI transition rate from m to g which 
will eventually be obtained from the present data anal- 
ysis. A variety of collision processes can lead to the 
quenching of m with simultaneous population increase 
of g. Among these processes are collisions of the stored 
ions with residual gas particles and quenching by external 
electromagnetic fields. In the above rate equations these 
processes are accounted for globally by the rate Aq. Ra- 
diative or dielectronic capture of a free electron by an 
initially metastable ions followed by Auger decay to the 
ground state, i.e., superelastic scattering (rate ^se) also 
leads to population transfer from m to g. Collisions may 
also result in excitation from g to m (rate A^) and in 
loss of ground-state and metastable ions from the stor- 
age ring (rates Agi and Ami, respectively). Finally, the 
rate Aur accounts for the depopulation of the metastable 
state by dielectronic recombination. At the experimen- 
tal clcctron-ion collision energy this loss process is only 
available for metastable ions and not for ground-state 
ions. 

With the initial conditions iV„(0) = Ng{0) = 

Ng°^ and with the definitions 



dt 



dt 



-{Ar + Ag + AsE + Adr + Ami)N„ 



+A^Ng, 



(1) 



= {Ar +Ag+ AS^)N^ - {Ae + Agl)Ng, (2) 



A = -{Ar + Ag + AsE + Ae + AuR + Ami + Agi) (3) 



T = ^J{A- Agif - Ae{AnR + Ami - Agt), (4) 
the solutions of Eqs. 1 and 2 are [42] : 



Nm{t) = I cosh(Tt) + [AeM°) -{A -A,- Agi)N^^' 



Ng{t) = I N^°^ cosh(Tt) + \{Ar + Ag+ Ase)N^^ + {A-Ae- A^O^i"' 



'^^]eM-At), (5) 
^^}exp(-^t). (6) 



Rewriting the hyperbolic functions as sums of two ex- 
ponentials reveals that the time evolution of the number 
of stored ions for both states is determined by the two 
time constants 



l/r„ 



A- 



T, 
T. 



(7) 
(8) 



which directly reflect the effective storage lifetimes of the 
metastable and ground-state ions, respectively, as sug- 
gested by the indices. As seen earlier [41], rather simple 
expressions for these time constants, are obtained in the 
limit = 0. In this case, T = A — Agi (Eq. 4), and Eqs. 
7 and 8 yield 



l/r™ 



Ar - 

Agi 



An + AsE + ^DR + A, 



ml 



(9) 
(10) 



As discussed in Sec. Ill C, this approximation is appro- 
priate in our case. Moreover, it will be shown that further 
relations between the parameters in these equations can 
be expected to be approximately fulfilled, so that the ra- 
diative lifetime Ar can be found once the decay times Tm 
and Tg have been experimentally determined. 

In the experiment, the measured count rate of re- 
combined ions comprises contributions from both the 
metastable and the ground state. Accordingly, it can 
be represented as the sum 



m 



-Pr, 



dN„,(t) 
' dt 



-Pg- 



d Ngjt) 
dt 



(11) 



where Pm (Pg) is the probability that a metastable 
(ground state) ion which is lost from the storage ring will 
then be collected by the recombination detector. This 
probability is effectively the rate of capturing an electron 
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either from clcctron-ion recombination in the Cooler or 
from the residual gas in the straight section of the stor- 
age ring leading up to the detector relative to the rate of 
all loss processes along the entire ring. 

Since the hyperbolic functions in Eqs. 5 and 6 can be 
expressed as the sums of two exponentials, the measured 
count rate (Eq. 11) can also be written as 



R{t) = Cme-^l^rn + Cge'^/^g 
where the coefficients 



Pm 

2t„ 



+ 



Pg 
It 



(l_^)Ar(o)_ -4- + -4. + -4sE ^(o) 



Co = 



2t„ 



(12) 



(13) 



(14) 



+ 



Pg_ 

2Ta 



with 



1, 



(15) 



are obtained by inserting Eqs. 5 and 6 into Eq. 11 and 
comparing the resulting expression with Eq. 12. For the 

lifetime determination described below, the coefficients 
Cm and Cg, obtained from fitting Eq. 12 to the data, will 
not be further considered. 



B. Fits to mesisured decay curves 

The results from the fits of Eq. 12 to the measured 
decay curves, i.e., the fit curves and the corresponding 
residuals, are included in Fig. 4. The best-fit values of 
the parameters r^, Tg, Cm and Cg are given in columns 
1-7 of Tab. IV along with their standard deviations. A 
range of initial data points were excluded from the fits in 
order to suppress the time intervals when the equilibrium 
conditions of the electron-cooled ion beam were not yet 
reached. The evolution of the best fit parameters when 
increasingly more early-data points are excluded from the 
fits is shown in Figure 6. For data set A, stable conditions 
were reached only after about 7 s (8 s for data sets B and 
C). Hence, data from the first 7 s (8 s) after ion injection 
into the TSR were disregarded in the fits to the three 
data sets, respectively. 

The error bars in Fig. 6 correspond to one standard 
deviation and were obtained from the nonlinear least- 
squares fitting procedure. These errors neglect mutual 
correlations between the fit parameters. The correlations 
can be evaluated from the covariance matrix of the fit [43] 
and be visualized as ellipses in two-dimensional subspaces 
of the fit-parameter space (Fig. 7a). In such plots the 
uncertainties of the individual fit parameters correspond 



34.0 

to 

32.0 - 



A=32 







i.O - 




7.0 



122.8 



122.7 



A=33 








6 8 10 12 14 4 6 

Start time of fit (s) 
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FIG. 6. Dependence of the fit parameters Tm (full symbols) 
and Tg (open symbols) on the time of the first data point 
that is used in the fit of Eq. 12 to the measured decay curves 
of data set A (Fig. 4a and 4b). The arrows mark the data 
that were used in the calculation of thfi- (Note that all values 
shown here scatter by less than the one-sigma error bars since 
they are statistically not independent, see text.) 
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FIG. 7. Correlations between the fit parameters Tm and Tg for 
data set A from the covariance matrices of the least squares 
fits to the experimental resonance-decay curve (Figs. 4a and 
4b). The ellipses denote the one sigma confidence limits for 
both Tm and Tg jointly. The dash-dotted horizontal and ver- 
tical lines represent the projections of the ellipses onto the 
parameter axes. 



to the projections of the ellipses onto the parameter axes. 
The resulting uncertainties are the ones that are given in 
Tab. IV. Since they stem from confidence regions for 
two parameters jointly, they are larger (in the present 
case by 50%) than the error bars in Fig. 6 that represent 
one-dimensional one-sigma confidence intervals. 
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TABLE IV. Isotope-dependent results for the decay times Tm,g and relative weights Cm,g obtained from the fits of Eq. 12 to the 
experimental decay curves displayed in Fig. 4. The standard deviations for these quantities were obtained from joint confidence 
intervals for Tm and Tg and for Cm and Cg, respectively (see text). Also given are the mass number A and the reduced values 
of the fits as well as the derived electron collision rates ^ec (Eqs. 18 and 19) and HFI lifetimes thfi (Eq. 21). 



data set 


A 




Tm. (s) 


Tg (s) 




Cm (S ^) 


Cg (s-^) 


^EC (s ^) 


Thfi (s) 


A 


32 
33 


1.010 
1.013 


33.0 ± 1.5 
8.0 ± 0.3 


136.3 ± 
122.7 ± 


0.4 
0.1 


451 ± 14 
722 ± 34 


3890 ± 17 
6565 ± 3 


0.023 ± 0.001 
0.022 ± 0.001 


10.5 ± 0.5 


B 


32 

33 


1.029 

1.026 


24.8 ± 2.0 
7.4 ± 0.6 


89.8 ± 

83.4 ± 


0.2 

0.1 


2226 ± 139 

1202 ± 130 


49024 ± 181 
14363 ± 10 


0.029 ± 0.003 

0.028 ± 0.003 


10.5 ± 1.2 


C 


32 
33 


1.042 
0.982 


30 ± 2 
7.8 ± 0.4 


148 ± 
177 ± 


6 
1 


62 ± 2 
100 ± 1 


63 ± 3 
84 ± 1 


0.027 ± 0.002 
0.025 ± 0.002 


10.2 ± 0.7 



C. Extraction of the hyperfine-induced decay rate 



With Aq = 0, Eq. 9 reads 



In general, the set of equations 12, 13, and 14 cannot 
be used to unambiguously derive all quantities of inter- 
est from the fits to the measured decay curves. However, 
the interpretation of the fitted decay constants in terms 
of atomic transition rates is facilitated if the effect of 
excitation of ground-state ions to the metastable state 
(rate A^) can be neglected, so that Eqs. 9 and 10 can be 
applied. CoUisional excitation has been investigated in 
previous experimental measurements of metastable state 
lifetimes, e.g., with stored C''+ [41], Xe+ [44], and Li+ 
[45] ions. In all these cases coUisional processes were 
found to have a negligibly small influence on the popula- 
tion of the metastable states under investigation, mainly 
because of the low residual particle density in the ul- 
trahigh vacuum of the storage ring. Since the residual 
gas pressure and the coUisional excitation cross sections 
of the present experiment and of the above mentioned 
studies are of the same order of magnitude, coUisional 
excitation processes are assumed to be negligible within 
the experimental uncertainty, and, accordingly, Eqs. 9 
and 10 are considered to hold. 

Deexcitation processes in the residual gas, denoted by 
Aq, are neglected here based on the same argument used 
above for the coUisional excitation processes {Ae) and 
we set Aq = 0. Since the rate Aq also accounts for 
the quenching of the metastable level by external mag- 
netic and electric fields these effects are implicity also 
assumed to be negligible. This assumption is justified 
in view of theoretical investigations of Stark quenching 
of the 2s 2p ^Pq level in Be-like ions [46] and of mag- 
netic quenching 2s 2p ^Pq level in the Be-like Ti^^+ [28]. 
Accordingly, quenching by electric fields [46] becomes no- 
ticeable only at field strengths which are orders of mag- 
nitude larger than the motional electric fields in the TSR 
dipole magnets. Likewise, in Ref. [28] it was shown that 
magnetic fields with strengths that are typically encoun- 
tered in the TSR do not significantly infiuence the life- 
time of the Ti^*''"(2s2p ^Pq) level. The present measure- 
ments which were performed at different settings of the 
TSR dipole magnets confirm this theoretical finding (see 
discussion below). 



=Ar+Aml + AuR + AsE- 



(16) 



For finding Ar the rates Ami , Aug and ^se still need to 
be specified. In order to provide a value for the rate Ami 
we make use of the plausible assumption that the coUi- 
sional rates of the ground and the metastable state are 
essentially equal, i.e., that A^i = Agi. This assumption 
is based on the fact that the ionization energies of the 
^So ground state and the ■^Pg metastable state, 652.2 eV 
and 627.5 eV, respectively [2], differ by less than 4%. It 
does not lead to any serious consequences. When, e.g., 
"^mi ~ ^^gf ^ assumed, the final result for the HFI 
transition rate changes by less than 1.5%. This change 
is smaller than the 4% experimental uncertainty (see be- 
low). 

Inserting Ami = Agi and Agi = I/tq (Eq. 10) into 
Eq. 16 yields 



1 

To 



A-EC 



(17) 



where the electron collision (EC) rate [47] Aec = ^dr + 
AsE accounts for dielectronic recombination of and su- 
perelastic scattering from initially metastable ions. Ex- 
ploiting the fact that there is no HFI transition for the 
A = 2,2 isotope, i.e., A^"^"* = 0, Eq. 17 yields the EC rate 
for this isotope: 



^EC 



1 



.^(32) 
/ m 



J32) • 



(18) 



The electron collision rate can be expressed as A^q = 

n^^^'uEC where the rate coefficient asc is the same for 
the metastable ions of both isotopes. Consequently, the 
EC rate for the A — ii isotope is 



,(33) 



J (33) _ ri_ 

^EC ~ (32) "^EC 



(32) 



,(33) 



1 



1 



,(32) I ^(32) 



J32) 



(19) 



In the previous experiment with Be-like Ti""^*"*" [27] the 
density ratio was very close to 1 and therefore omitted 
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from the equations. In the present measurements the 
ratios were 0.938. 0.948, and 0.932 for data sets A, B, 
and C, respectively (Tab. III). As can be seen from the 
next to last column of Tab. IV the rates Aec for the three 
experimental data sets mutually agree within 5%. 

Using now Eq. 17 for the ^ = 33 isotope gives the 
2s 2p ^Pq — >• 2s^ transition rate in the laboratory 
frame: 



^(33) ^ 



,(3.-5) 



1 



.(33) 



.(33) 



n, 



(32) 1 (32) 
\ ' m 



.(32) 



(20) 



Finally, the HFI transition rate is obtained from a trans- 
formation into the reference frame of the moving ions, 



(33) 

i.e., by multiplication of Ai- with the Lorentz factor 
(33) = [1 _ (^(33)/^)2]-i/2 denotes the speed of light): 

1 
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iHFI 



thfi 



(21) 



■7 



(33) 



,(33) 



J33) 



J33) 



(32) 1 ^(32) 



J32) 



Results for the hyjxafine induced lifetime as derived 
from the three experimental data sets A, B, and C, are 
given in the last column of Tab. IV. This table also lists 
the best fit values for the decay time constants which were 
used in Eq. 21 together with the experimental parame- 
ters from Tab. III. Although the three data sets were 
obtained under significantly different experimental con- 
ditions, the three experimental values for thfi agree with 
each other within the experimental uncertainties. This 
documents the robustness of the present experimental 
approach. The weighted mean (with the inverse squared 
errors as weights) of the three experimental HFI lifetimes 
is Thfi = 10.4 ± 0.4 s. 



IV. DISCUSSION 

During the measurement of data set C the magnetic 
field in the TSR bending magnets was a factor of two 
stronger as compared to the measurements of data sets 
A and B (Tab. III). At the same time, thfi from data 
set C is smaller by 3% than thfi from data sets A and B 
(Tab. IV). This may be an indication of a slight quench- 
ing of the 2s 2p ^Pq state in the magnetic fields of the 



TABLE V. Comparison of theoretical values for the hyperfine 
induced 2s2p ^Pq —¥ 2s^ ^So transition rate in Be-like S^^"*" 
with the present experimental result. 



Year 


^HFI 

Experiment 


Theory 


Deviation from 
experiment 


Ref. 


1993 




0.03611 


-62.4% 


[6] 


2008 




0.09315 


- 3.0% 


[8] 


2009 




0.09355 


- 2.5% 


[9] 


2011 


0.096(4) 






present 



storage ring magnets (sec also [48]) which cover 13% of 
the closed orbit of the stored ions (Fig. 2). This effect 
is expected to increase linearly with increasing magnetic 
field strength. Since the 3% change of thfi is smaller 
than the 4% experimental uncertainty, it must, how- 
ever, be concluded that there is no significant influence of 
the TSR magnetic fields on the measured HFI transition 
rates at the present level of experimental accuracy. This 
conclusion is in accord with recent theoretical results for 
isoelectronic Ti^^+ [28]. 

In Tab. V the present experimental result for the 
HFI transition rate ^hfi = I/^hfi = 1/10.4(4) s"^ = 
0.096(4) s~^ is compared with the available theoretical 
results [6, 8, 9]. The early theoretical calculation by Mar- 
ques et al. [6] yielded a value that is about 60% lower 
than the experimentally determined rate. This is in con- 
trast to the later calculations by Cheng et al. [8] and 
Andersson et al. [9] which both yield results in very good 
agreement with the present experimental value. As al- 
ready discussed (Sec. I), the main difference between the 
early and later calculations is the treatment of electron- 
electron correlation effects. In particular. Marques et 
al. [6] disregarded the mixing of the 2s 2p ^Pi level with 
the 2s 2p level. This mixing was included in the later 
calculations of HFI transition rates in Be- like ions [7-10]. 
The comparison of the available theoretical results with 
the present experimental result for 33gi2-i- strongly sug- 
gests that the additional computational effort of the later 
calculations is indeed required for arriving at more accu- 
rate HFI transition rates. 

Generally, the relative importance of correlation effects 
decreases with charge state. In the present context, this 
can be seen from the differences between the early and 
the later theoretical HFI transition rates. For the lower 
charged 33gi2+ result of Marques et al. [6] devi- 

ates by 61% from the result of Cheng et al. [8] (Tab. V). 
For the more highly charged 47xii8-H 

ion this difference 

is only 47% (Tab. I). In view of this fact, the ~20% 
disagreement of the experimental result for 47-^^184- .^^i^h 
the latest theoretical calculation [8, 9], which agree well 
with the 33gi2-i- experimental value, is surprising. In the 
attempt to find an explanation for this disturbing dis- 
crepancy, the influence of external magnetic flelds on the 
lifetime of the '^'^Ti^^+{2s2p ^Pq) has been investigated 
by detailed theoretical calculations [28]. From this study 
it was concluded that the magnetic fields in the the TSR 
bending magnets do not influence the measured HFI life- 
times significantly. This is in agreement with the present 
experimental findings for 33gi2-i- jQjjg_ 



V. CONCLUSIONS 

We have measured the hyperfine induced 2s 2p ^Pq 
2s^ ^So transition rate in Be-like 33gi2-i- -^j^g ^gjj^g jgQ. 
topically pure ion beams in a heavy-ion storage ring. The 
decay of the number of stored ions was monitored by 
employing electron-ion recombination. In order to selec- 
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tivcly enhance the signal from the 2s 2p ^Pq level, the 
electron-ion collision energy was tuned to a dielectronic 
recombination resonance associated with this level. 

In our present experimental work, different data sets 
were obtained under substantially different conditions. 
Nevertheless, the individual HFI transition rates from 
the different data sets agree with each other within the 
experimental uncertainties. This demonstrates the ro- 
bustness of the present experimental approach. 

Our experimental value of 0.096(4) s~^ agrees with the 
latest theoretical calculations [8, 9]. No significant influ- 
ence of the magnetic field in the storage-ring bending 
magnets on our result has been found. The origin of the 
discrepancy between our previously measured HFI rate 



for 4^Tii^ [27] and the latest calculations [8-10] still re- 
mains unclear. 
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